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The recently reported precise experimental determination of the dipole polarizability of the H^ 
molecular ion ground state [P.L. Jacobson, R.A. Komara, W.G. Sturrus, and S.R. Lundeen, Phys. 
Rev. A 62, 012509 (2000)] reveals a discrepancy between theory and experiment of about 0.0007ag, 
which has been attributed to relativistic and QED effects. In present work we analyze an influence 
of the relativistic effects on the scalar dipole polarizability of an isolated H^~ molecular ion. Our 
conclusion is that it accounts for only 1/5 of the measured discrepancy. 



I. INTRODUCTION 
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Recent measurements fl|J2j of the scalar electric dipole polarizability of Hj molecular ion through the study of 
H2 molecule states with one Rydberg electron stimulated the introduction of methods j|- || which are able to 
accurately describe wave functions of molecular ions with two heavy nuclei beyond the adiabatic approximation. 
The accuracy for the dipole polarizability constant (~ 10~ 7 aQ) reached in the last work M in its turn become a 
challenge to experiment. The new experimental work B substantially increases the accuracy of measurements and 
reveals a discrepancy of about 0.0007a 3 ) between theory and experiment, which can not be accounted for within 
purely nonrelativistic approximation. In present work we consider relativistic corrections of order a 2 to the dipole 
polarizability of the ground state of an isolated HJ" molecular ion. 
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II. THEORY 



The nonrelativistic Hamiltonian of the hydrogen molecular ion HJ" is 
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We adopt atomic units (e = % = m = 1) throughout this paper. The interaction with an external electric field (details 
of the nonrelativistic treatment of the problem can be found in previous papers [J3|- fel) is expressed by 



Vn = £ n • d, 



(2) 



where 
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\2M + m 2M + m 
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is the electric dipole moment of the three particles with respect to the center of mass of the system. Without loss of 
generality we assume that n ■ d = \xz c . 

The Breit a 2 correction to the nonrelativistic Hamiltonian is described by an operator 



Vb=o?\- V 



U [S(r 1 )+S(r 2 )} + 



8m 3 8m 2 
Then the total Hamiltonian reads, 

H = H + V B + V P . 
Let us define the ground state nonrelativistic wave function as follows 
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{H - £ )*o = °- 
In the nonrelativistic case the change of energy due to polarizability of molecular ion is expressed by 
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and 



E^ = (* \V p (Eo - H )- l V p \* ) = E 2 /j, 2 {$ \z c (E - H )- X z c \* t 



a° s = -2n 2 (V \z c (E - H )- l z c \* ). 
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Let us introduce H\ — H + Vb , then the scalar dipole polarizability a s with account of relativistic corrections can 
be rewritten in a form (we assume that Vb ~ o?H§ and ^ B = "Jo + ^> B ) 



(8) 



a£= -2 f i 2 (^ B \z c {E 1 - H^z^Z) 

= -2fj, 2 {^ B \z c [(E - Ho)- 1 + {Eo - Ho)- 1 (V B ~ (Vb)) {Eq - Ho)' 1 + . . .] z c \V B ) 
= -2 M 2 (* |* C (£ - Ho^z^o) - 2 f i 2 (^o\z c (E - H )- 1 V B (E Q - flb) _1 «c|*o) 
-2 f i 2 ^ B \z c (E - F )-^c|*o) + {*o\zc(Eo - H )"^c|* s )) , 

and $> B = (E Q — H )~ 1 Vb\^o)- Thus relativistic correction to the scalar dipole polarizability a s is reduced to 
evaluation of the following matrix elements 

A« s =-2fi 2 {^ \z c (E Q - Ho)- 1 {V B - (Vb)) (E - flb) _1 *c|*o) 

-2^ 2 {^o\V B (Eo - Ho)- x z c (Eo - H )- l z o \^ ) (9) 

-2fi 2 (y \z c (Eo - H )- l z c (E - J?o)- 1 V B |*o). 

At this point we can note that the spin-orbit term does not contribute to a s since the magnetic dipole operator has 
selection rules m! = m ± 1. 

III. VARIATIONAL NONRELATIVISTIC WAVE FUNCTION 



Variational wave function describing the ground state of the hydrogen molecular ion is taken in a form 
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| Ci cos{viR 12 ) + A sin(viRi 2 ) 
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xe -* ,i "*" """ + (ri <-> r 2 ). 
Here a%, Pi, 7j, and Vi are parameters generated in a quasirandom manner, 



1., . _ 



{(A 2 -A 1 ) + A x ], 



where \_x\ designates the fractional part of x, p a and q a are some prime numbers, the end points of an interval [A\ , A 2 ] 
are real variational parameters. Parameters Pi, ji, and Vi are obtained in a similar way. Details of the method and 
discussion of various aspects of its application can be found in our previous papers |7]|| ■ 
The perturbed function \&i = \i{E§ — Hq) r c ^o is expanded in the similar way 



*i =J^ ri Ci cos(uiR 12 ) + A sin(^i?i 2 ) 
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Technically evaluation of the matrix elements in Eq. (|9|) can proceed in the following way. For the matrix element 
in the first line we need to solve one linear equation, (Eq — -ffo)^i = z c^o, and then average operator (Vb — (Vb) 
over $i. To get rid of singularities in the solutions of implicit equations of lines 2 and 3 of Eq. ([))) one can solve a 
sequence of equations from the right to the left in the second line and in the reverse order for the third line. In the 
latter two cases solution as well as the rhs of the last equation should be projected onto subspace orthogonal to \&o- 



IV. RESULTS AND CONCLUSIONS 

In Tabic I results of numerical calculations are presented. For the zeroth-order approximation a wave function with 
a basis set of N = 800 has been used that yields the nonrelativistic energy 

E = -0.59713906312340(1), (12) 

which is in a good agreement with our previous accurate result [g|. Here m p — 1836.152701m e is adopted. As seen 
from the Table convergence for the relativistic correction is slower due to singular operators encountered in the matrix 
elements. Nevertheless we can conclude from this Table that the resulting value is 

Aa s = -0.00015214(1). (13) 

Combining this result with the nonrelativistic value from the Table I one obtains that the static electric dipole 
polarizability of H.j~ molecular ion ground state with relativistic a 2 corrections is to be 

a s = 3.1685737(1). (14) 

We see that thus obtained value for the dipole polarizability does not fully account for present disagreement between 
theory and experiment (comparison of our results with results of previous calculations and experiment are presented 
in Table II). Our consideration does not include leading order QED corrections but usually they are one order of 
magnitude smaller than relativistic corrections and have a different sign. Thus they could not cover the rest 4/5 of 
the discrepancy. 

On the other hand the experimental value for the dipole polarizability has been deduced from the effective model 
Hamiltonian Jll[ ] which is a fully nonrelativistic Hamiltonian and it does not include the retardation Casimir-Polder 
effect |12|| f or the Rydberg electron. On the importance of this phenomena has been pointed out in a paper of Babb and 
Spruch|13[| . So, our conclusion is that the Casimir-Polder potential has to be included into the effective Hamiltonian 
to meet the requirements of the present level of experimental accuracy. That will enable to deduce scalar electric 
dipole polarizability in a more reliable way. 

This work was supported by the National Science Foundation through a grant for the Institute for Theoretical 
Atomic and Molecular Physics at Harvard University and Smithsonian Astrophysical Observatory, which is gratefully 
acknowledged. 
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«2 Aa 3 

AT = 400 3.1685962 -1.52065 [-4] 

TV = 600 3.1687252 -1.52140[-4] 

TV = 800 3.1687258 -1.52137[-4] 

TABLE I. Dipole polarizability of H^(0, 0). Convergence of the numerical results with the size of basis set. 



Experiment [| 3.16796(15) 

Shertzer and Greene || 3.1682(4) 

Bhatia and Drachman @ 3.1680 

Moss |] 3.16850 

Taylor, Dalgarno, Babb § 3.1687256(1) 
present work 

nonrelativistic 3.1687258 

with a 2 corrections 3.1685737 

TABLE II. Dipole polarizability of Hj(0, 0). Comparison with other calculations and experiment. 



